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Delocalization of Nonlinear Optical Responses in Plasmonic Nanoantennas 
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Remote excitation and emission of two-photon luminescence and second-harmonic generation 
are observed in micrometer long gold rod optical antennas upon local illumination with a tightly 
focused near-infrared femtosecond laser beam. We show that the nonlinear radiations can be emitted 
from the entire antenna and the measured far-held angular patterns bear the information regarding 
the nature and origins of the respective nonlinear processes. We demonstrate that the nonlinear 
responses are transported by the propagating surface plasmon at excitation frequency, enabling 
thereby polariton-mediated tailoring and design of nonlinear responses. 

PACS numbers: 78.67.-n, 78.67.Qa, 78.67.Uh, 78.60.Lc, 42.65.Ky 


Optical antennas are pervasive devices to control 
spatial distribution of light on sub-diffraction length 
scales [1-3]. Concurrently it is realized that held en¬ 
hancing properties of underlying surface plasmon (SP) 
resonances may foster much needed nonlinear behaviors 
to improve nanoscale light management [4-6]. A nonlin¬ 
ear optical antenna combines the functionalities of linear 
devices (extreme light concentration, tailoring of spatial 
and phase distributions, directivity of emission, etc) with 
the benefits of nonlinear optical effects, such as frequency 
conversion [7], mixing [8], ultrafast switching, modula¬ 
tion [9, 10] and self-action effects [11, 12], to name just 
a few. Nonlinear responses are notably intricate and 
a comprehensive theory of nonlinear plasmonics is still 
underway. To date the most extensively elaborated are 
harmonic generation, four-wave mixing and multi-photon 
luminescence at the single optical nanoantenna level [13- 
22]. It has been repeatedly demonstrated that nonlinear 
spectral and intensity responses in plasmonic antennas 
are largely determined by the localized SP resonances 
at the frequencies of the driving optical fields [5, 23-25]. 
In spatially extended plasmonic objects, point-and-probe 
nonlinear scanning microscopy revealed the importance 
of the supported SP modal landscape [26-30], suggesting 
that the nonlinear responses may bear the signatures of 
the SP mode spatial extension. It is this effect that we 
address in this letter. 

To this aim, we discuss two nonlinear processes - two- 
photon luminescence (TPL) and second-harmonic gener¬ 
ation (SHG) - from gold rod optical antennas upon lo¬ 
cal illumination with a tightly focused femtosecond near- 
infrared laser beam. We show that in this type of struc¬ 
tures, nonlinear confocal TPL and SHG mappings are in¬ 
effective when it comes to discerning differences between 
these two processes. The variations between incoher¬ 
ent TPL and coherent SHG are unambiguously revealed 
in spectrally filtered angular distributions measured in 
Fourier and image planes. Importantly, we demonstrate 
that nonlinear conversions of the incident electromag¬ 
netic energy are not restricted to the excitation area but 
are spatially delocalized along the entire structure. We 


argue that nonlinear optical transport is mediated by a 
propagating SP at the excitation frequency despite the 
associated high losses. We substantiate this hypothesis 
by modelling far-held SHG signatures as originating from 
a line of non-interacting dipolar sources oscillating at the 
SHG frequency, whose amplitudes and phases are deter¬ 
mined by the damped SP polariton at the fundamental 
frequency, developing in the one-dimensional cavity-like 
antenna. Our findings demonstrate new degrees of free¬ 
dom for design of SP mediated coherent and incoherent 
nonlinear processes. 

Optical gap antennas are fabricated by electron-beam 
lithography and lift-off technique on a glass substrate. 
Each antenna consists of two identical nanowires sepa¬ 
rated by a gap. The dimensions of individual wires are 
110 nm in width and 50 nm in height. The length L of 
each arm and the gap separation g are systematically var¬ 
ied from 300 nm< L <4000 nm, and from 0 nm< g <150 
nm, with a minimum gap size of ca. 20 nm, as mea¬ 
sured with scanning electron microscope (SEM). Optical 
excitation and collection are performed using an inverted 
microscope. A 180 fs Titanium: Sapphire laser tuned at a 
fundamental wavelength of Ao =810 nm is focused on the 
antennas in a diffraction-limited spot by a high numerical 
aperture objective (oil immersion, NA=1.49). The aver¬ 
age laser power at the sample is 3.5 mW. The incident 
beam is linearly polarized along the antenna. Nonlin¬ 
ear signals are collected by the same objective followed 
by a dichroic beam splitter, which separates the useful 
spectral range (375-700 nm) from the back-scattered fun¬ 
damental beam. Simultaneous TPL and SHG confocal 
maps are collected by two avalanche photodiodes in the 
absence of a spatial filter, allowing detection of signal 
emitted from the entire structure. A 10 nm narrow band¬ 
pass filter centered at 405 nm is used for SHG detection. 
Fourier and direct planes imaging are recorded by sep¬ 
arate cameras and relay lenses appropriately placed in 
their respective conjugate planes [31]. 

TPL is an incoherent nonlinear optical process defined 
by the material’s electronic band structure [32-34] and 
underlying plasmonic modes [15, 25, 29, 35]. In con- 
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FIG. 1. (a)-(c) TPL, SHG confocal maps and SEM micro¬ 
graphs of three gold rod antennas. The arm lengths of the 
coupled antennas are 650 nm, the gaps are g=130 and 50 nm. 
(d) and (e) Angular distributions of TPL and SHG signal 
from a single nanowire (ca.1500 nm long, the bottom antenna 
in (a)-(c)). Scale bars are 200 nm. The laser is focused at 
the left extremity. The images are normalized in units of 
kem = 27r/A e m, where A em is the emission wavelength of the 
nonlinear process. The inter-fringe distance Ak is denoted 
with the white double-headed arrow, (f) Linear dependence 
of 2iv/Ak x on nanowire length L (solid line is the linear least 
square fit) and 2n/Ak x = L (dotted line). 


junction with confocal mapping technique, it probes lo¬ 
cal electric field intensity in plasmonic objects [26, 36]. 
Similarly, SHG point-and-probe mapping provides some 
additional information about material due to the intrin¬ 
sic dependency of SHG on structural symmetries [27, 28, 
37, 38]. Having applied both methods to the lithographic 
rod antennas, we find that, apart from the difference in 
the nonlinear yield (typically TPL/SHG ratio is at least 
10 3 ), confocal nonlinear maps are practically identical as 
illustrated in Fig. 1(a) and (b). This was true for all 
studied optical gap antennas, regardless of the length or 
the gap size. Because of this striking similarity, we rule 
out any significant structural symmetry-dependent SHG 
component, which could render SHG maps to be some¬ 
what different from the TPL ones [5]. It follows that the 
spatial variations of both SHG and TPL signals must be 
accounted for by the same sensitivity to the local plas¬ 
monic modal distribution [29, 39]. 

While nonlinear confocal mapping fails at distinguish¬ 
ing between these fundamentally different nonlinear pro¬ 
cesses, Fourier plane imaging unveils the information hid¬ 
den in the pixels of the nonlinear confocal maps. In 
Fig. 1(d) and (e), we compare Fourier plane images rep¬ 
resenting the projected angular distributions of the TPL 
and SHG emissions from a single rod nanoantenna [bot¬ 
tom nanowire in Fig. l(a)-(c)]. When the laser is fo¬ 
cused on the left extremity, the TPL Fourier distribution 
[Fig. 1(d)] features a pattern with two maxima aligned 
along the antennas x axis coinciding with the excitation 
polarization direction at variance with a single x-oriented 
dipole [40]. Vastly different is the SHG angular image 


shown in Fig. 1(e) displaying distinct interference fringes 
reflecting the coherent nature of SHG. Excitation at the 
left extremity of the antennas systematically produces 
fringe patterns with an intensity increasing towards pos¬ 
itive k x /k em - Such fringe pattern was already observed 
in the linear regime for nanowires decorated with emit¬ 
ters [40]. SHG Fourier images show a strict dependence 
on the rod length. Figure 1(f) shows the reciprocal de¬ 
pendence of the fringe period A k x on the antenna length. 
As TPL is intrinsically an incoherent process, we do not 
observe interferences in Fourier plane even when a nar¬ 
row portion of the broad TPL spectrum was spectrally- 
selected. Neither do the TPL Fourier images depend 
on the antenna’s length and patterns similar to that in 
Fig. 1(d) are persistently obtained. 

Fourier plane imaging of SHG response of coupled op¬ 
tical gap antenna shows sensitivity to the gap size. Fig¬ 
ure 2(a)-(c) illustrate a set of Fourier planes obtained 
from three antennas having identical arm lengths (ca. 
830 nm) and gaps of 140 nm, 40 nm and at contact, cor¬ 
respondingly. The excitation is located at their respective 
right extremities. For maximally decoupled arms [<?=140 
nm, Fig. 2(a)], we observe four wide fringes, which, as 
the gap size decreases, start to split [Fig. 2(b)], forming 
a pattern of a different symmetry. The splitting becomes 
more pronounced as the gap decreases to its minimum. 
A Fourier pattern corresponding to the case of touching 
arms [Fig. 2(c)] contains a double number of fringes as 
compared to the case of maximally separated arms. We 
attribute such splitting/pairing of the fringes to the onset 
of coupling between the antenna’s arms. 

The sensitivity of SHG Fourier planes to the antenna’s 
length and the gap size indicates that the nonlinear re¬ 
sponse is not simply generated locally by the focused 
laser beam. At this point, we note that the excita¬ 
tion of the nanorod extremity by a focused laser beam 
is a prevalent technique to launch a propagating SP in 
the structure [41-43]. In order to account for the SHG 
Fourier planes’ length dependence, in the following we 
assume that the propagating SP is creating an enhanced 
electric field along the antenna strong enough to allow 
remote nonlinear optical interactions. We construct a 
model of the delocalized SHG emission from a single gold 
antenna of length L by considering propagation of the 
SP in one-dimensional cavity (ODC) of length L s i m [44]. 
Upon point dipole excitation at the cavity’s left extrem¬ 
ity (x = —L s i m /2), which we equate here with the fo¬ 
cused laser beam excitation, the plasmon associated elec¬ 
tric field along the antenna Ef (x) can be written as [44]: 


E£ (x ) = E 0 e ikspL ‘ im/2 (r+l)- 


o^SI 


_ rppi'ks'p (L s: 


')) 


1 — r 2 e 


2 zjz2fcsp L 

s 


(1) 


where Eq is the excitation field amplitude, r is the reflec¬ 
tion coefficient. The values of complex SP wave vector 
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FIG. 2. (a)-(c) Evolution of Fourier plane patterns with the 
gap size: 140 nm, 40 nm, 0 nm. Antennas arm length is ca. 
830 nm, excitation at the right-hand extremity. Simulated 
Fourier planes: (d) L s im=500 nm, (e) far-held beating be¬ 
tween L s i m =500 nm and 1400 nm, (f) L s i m = 1400 nm. The 
effective refractive index at to is n e fr =1.73, the propagation 
length is L sp = 1/2 k" p = 2100 nm. 


fc sp = fc'p + ik' s ' p are extracted from finite element simula¬ 
tions of an infinitely long nanowire with the cross section 
of the nanowires used in the experiment. The nonlin¬ 
ear emission is phenomenologically modelled by a large 
number of identical non-interacting effective dipoles p 2u , 
oscillating at the SHG angular frequency. The dipoles 
are aligned and oriented along the antenna x axis on the 
glass/air interface. The amplitudes and phases of the 
individual oscillators are position-dependent and deter¬ 
mined by the ODC scalar plasmon wave E“ p (x), accord¬ 
ing to p 2ul (x) = f3^[Es P (x)] 2 , where is a ID counter¬ 
part of the nonlinear polarizability of the individual ef¬ 
fective dipoles [45]. This treatment of SHG is inspired by 
the earlier works establishing the link between the dom¬ 
inant SHG origin with the local symmetry breaking at 
the surface [16, 17]. Note that in our case the interband 
structure of gold inhibits plasmon modes at the SHG fre¬ 
quency due to strong absorption and prevents tailoring 
of plasmonic modes at the harmonic frequency [46]. 

We employ the known expressions of the far-field Green 
functions to calculate the coherent sum of electric fields 
emitted by each SHG dipole source, collected by the ob¬ 
jective [47] . The results of simulations for the line of p 2ul 
dipoles with the total length of L s i m =1400 nm and 500 
nm are shown in Fig. 2(d) and (f). In the presence of a 
gap, we calculate the Fourier plane as a coherent beating 
between far-held signals emitted from two nanowires of 
-b s im=500 nm and £ s im=1400 nm. The result of the cal¬ 
culation is shown in Fig 2(e), which reproduces well the 
experimental fringe splitting when short (S) and long (L) 


lengths contribute to the final signal with equal weights 
ws = wl- We find a good match between experimental 
and simulated data when L sim is close to the effective 
nanowire length 27 t/A k x [Fig. 1(f)]. We think that the 
systematic difference between the physical length L and 
Lsim could be due to intrinsic oversimplification of the 
real electric held in the nanowires inherent to the OCD 
model. Among other factors is the omission of other pos¬ 
sible mechanisms of SHG. The latter will be particularly 
interesting to investigate in the light of the recent discus¬ 
sions on ponderomotive force, Kerr-like and heat induced 
nonlinearities in metals [48, 49]. 

To accurately ht the experimental SHG Fourier images, 
the coefficient r in Eq. 1 is set to zero. This implies that 
there is no contra-propagating held components in the 
cavity, which can be understood by large Ohmic losses 
and efficient end-face scattering [42]. Interestingly, our 
simulations predict that if even a small portion of the 
electromagnetic excitation could scatter back into the 
cavity, a drastic change would be observed in the SHG 
Fourier pattern, which is not the case here. In the ab¬ 
sence of a back-reflected SP, the ODC model prediction 
of the scalar held E£ p (x) coincides simply with a damped 
wave e* fes p( :c + I ' si ”/ 2 ), as can be easily seen from Eq. 1. In 
this case, the Fourier planes are merely Fourier trans¬ 
forms of the e 2 * fc sp( x + i sim/2) f urLC tion, and the observed 
fringes are so-called Gibbs oscillation, resulting from the 
hnite length of waveguide [50]. 

Figure 3 experimentally conhrms the assumptions of 
delocalized nonlinear responses tested in the simulation. 
Figure 3(a) shows an image plane micrograph of a single 
rod antenna excited locally by the laser beam focused at 
its left extremity. The image is recorded at the funda¬ 
mental wavelength Ao- The light detected at the distal 
end unambiguously indicates the excitation of a SP mode 
in the antenna [41]. No leaky plasmon mode is expected 
in our structures due to their small transverse cross sec¬ 
tion [42]. The corresponding spectrally hltered TPL and 
SHG image planes are shown in Fig. 3(b) and (c), re¬ 
spectively. Aside from the strong local spot at the laser 
position, it becomes evident that the TPL response is 
delocalized along the entire rod [Fig. 3(b)]. Local scat¬ 
tering at the rods defects and structural discontinuities 
are also readily observed. Similar point scattering is ob¬ 
served in the SHG filtered image of Fig. 3(c), together 
with characteristic interference patterns reminiscent of 
the fringes detected in Fourier planes. While these non¬ 
linear images are somewhat degraded by residual chro¬ 
matic aberrations from relay lenses, they unambiguously 
demonstrate the SP mode developing in the antenna at 
Ao carries enough energy to produce a distributed non¬ 
linear response during its propagation. To further con¬ 
firm this conclusion, we studied the delocalized responses 
for coupled gap antenna with g =50 nm and g=180 nm. 
The images recorded at Ao indicate a significant near-field 
coupling for the smallest gap as the plasmon excitation 
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FIG. 3. (a), (b) and (c) are colorized direct image plane mi¬ 
crographs emphasizing the excitation of propagating SP mode 
in a 1950 nm-long gold rod nanoantennas. Images are spec¬ 
trally filtered to record emission at the fundamental, the TPL, 
and SHG wavelengths, respectively. The nonlinear signals are 
delocalized along the entire length of the antenna, (d) to (f) 
and (g) to (i) are obtained from coupled gap antennas with 
g= 50 nm, and g=180 nm. Insets are SEM images of the cor¬ 
responding antennas. Excitation is at the extreme left end. 

is transferred from one arm to the other [Fig. 3(d)]. For 
the decoupled antenna, scattering at the gap strongly 
mitigates the transmission to the right arm [Fig. 3(g)]. 
Concomitantly, the spatial extent of the nonlinear pro¬ 
cesses are gap dependent. In Fig. 3(e)-(i), both TPL and 
SHG emission are observed from the second arm, in line 
with the fact that the plasmon at Ao is conveyed through 
the gap. For decoupled antennas, the intensity of the 
nonlinear emissions emitted along the second arm disap¬ 
pear because most of the energy carried by the plasmon 
is scattered at the gap [Fig. 3(h)-(i)]. 

Hence, the efficiency of energy transport through the 
gap upon excitation at the extremity can be indirectly 
monitored by (i) observing light scattering from the gap 
in linear regime [Fig. 3(d) and (g)], or by (ii) observ¬ 
ing the fringe patterns in Fourier planes in nonlinear 
regime (Fig. 2). These two aspects are brought together 
in Fig. 4. The relative scattering efficiency at Ao is mea¬ 
sured as a function of gap size (orange circles) from an¬ 
tennas shown in Fig. 3, where maximal coupling corre¬ 
sponds to the minimal scattering (g= 0). Then, we plot 
the quantity ws/(ws + wl) describing the mixing of the 
contributing antenna’s dimensions. The black circles and 
stars in Fig. 4 are the data point obtained from 1950 nm 
and 830 nm arm length, respectively [Fig. 2 and Fig. 3]. 
The quantities ws and wl are inferred by comparing 
and matching the respective experimental and theoret¬ 
ical Fourier planes. It becomes apparent that the black 
and orange markers form a single trend. This correla¬ 
tion is a demonstration of the SHG Fourier plane imag¬ 
ing potency to monitor near-field coupling in plasmonic 
systems, and, in conjunction with theoretical modelling, 



g(nm) 

FIG. 4. Correlation between relative linear scattering efficien¬ 
cies of SP from the gaps of the 1950 nm arm long antenna 
(orange circles, Fig. 3(a),(d) and (g)) and the extracted from 
the experimental and simulated Fourier plane images values 
of the mixing parameter ws/(ws + wl) (black circles). Black 
stars are the mixing parameters for the set of ca. 830 nm 
arm long antennas (Fig. 2). The insets are experimental and 
simulated Fourier planes, showing different degrees of fringe 
splitting for two selected gap-antennas. The orange curve is 
added for eye-guiding purposes. 


to constitute a tool for nonlinear optical far-field recon¬ 
struction of plasmonic near-fields. 

In summary, by implementing spectrally filtered 
Fourier and image plane measurements and a simple phe¬ 
nomenological model, we demonstrate that nonlinear re¬ 
sponses in gold rod optical antennas can have a signifi¬ 
cant spatially delocalized contribution when excited lo¬ 
cally by a focused laser beam. This result presents an 
alternative - delocalized - mechanism of nonlinear emis¬ 
sion excitation in plasmonic structures, which can be 
adopted to observe and explain a variety of parametric 
and non-parametric nonlinear processes. We anticipate 
that incorporation of the delocalization concept in the 
current nonlinear plasmonic discourse can incite novel 
approaches for coherent control and tailoring of nonlin- 
earities in single plasmonic cavities and nretamaterials 
via plasmon mediated (quasi) phase matching. 
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